Ground states of the Josephson vortex lattice in layered superconductors 
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We consider the ground state configurations of tiie Joseplison vortex iattice in iayered supercon- 
ductors. Due to commensurability effects witfi ttie layered structure, tiie iattice tias multiple con- 
figurations, both aligned with layers and rotated at finite angle. At low fields the lattice switches 
between these configurations via first-order phase transitions. These transitions become more fre- 
quent at smaller fields. With increasing magnetic field a dilute lattice transforms first into a sheared 
dense lattice. With further increase of field, the shear deformation smoothly vanishes at a second- 
order phase transition. 



A magnetic field applied parallel to the layers gen- 
erates a lattice of Josephson vortices with many un- 
usual static and dynamic properties. An important field 
scale is set by the anisotropy factor 7 and interlayer pe- 
riodicity s, Bcr = $o/(27r7s^) (in Bi2Sr2CaCu208-<5, 
Bcr ^ 0.3 — 0.5T). There are two very different regimes 
depending on the strength of the magnetic field. In the 
dilute lattice regime, <^ Bcr, the nonlinear cores 
of Josephson vortices are well separated and the distri- 
bution of currents and fields is very similar to that in 
continuous anisotropic superconductors. At high fields, 
Bx > Bcr, the dense lattice regime is realized where the 
cores of the Josephson vortices overlap. In this regime the 
Josephson vortices fill all layers homogeneouslji. This 
state is characterized by very weak modulations of the 
in-plane and Josephson currents. In this proceedings we 
consider the evolution of the ground state configuration 
with increasing magnetic field. 

As the centers of the Josephson vortices must be lo- 
cated between the layers, the layered structure plays 
crucial role in selection of the ground-state lattice con- 
figurations. The Josephson vortex lattice is commen- 
surate with the layered structure only at the discrete 
set of magnetic fields. At small magnetic fields, B^ < 
<I>o/(27r7s^), the Josephson vortex lattice can be de- 
scribed by anisotropic London model (see, e.g., Refs. 
013) • This model has a simple scaling property: in scaled 
coordinates, z = z/s, y = y/^s, the energy becomes 
isotropic in z-y plane, which means that the ground state 
in these coordinates corresponds to an ideal triangular 
lattice which is degenerate with respect to rotations in 
this plane. In real coordinates these rotations correspond 
to "elliptic rotations". As a consequence, the family of 
commensurate lattices includes lattices aligned with the 
layers, as well as misaligned ones. To make a full clas- 
sification of commensurate lattices we consider a general 
lattice shown in Fig. QJ^^. The lattice is characterized 
by three parameters: in-plane period a, distance between 
vortex rows in c direction b = Ns, and relative shift be- 
tween the neighboring vortex rows in c direction qa. The 
lattice shape is characterized by the two dimensionless 
parameters, q and ratio r = b/a. The lattice param- 
eters are related to the in-plane magnetic field, B^, as 
B^ = ^o/i^b). Structures aligned with the layers cor- 
respond to 5 = 1/2. As the replacement q ^ 1 — q 
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FIG. 1: (a) General Josephson vortex lattice and its parame- 
ters, (b) Orientation of layered structure with respect to ideal 
lattice (in scaled coordinates). Layered structure fits the ideal 
lattice only if it is oriented along one of the crystallographic 
directions, which is characterized by two numbers (n, m), cor- 
responding to expansion of the direction vector over the two 
basic lattice vectors, ei and 62. Several possible directions 
are shown with the corresponding indices (n,m). The lower 
part of the figure illustrates the lattice parameters, a, 6, and 
q for one possible orientation. 



corresponds to the mirror reflection with respect to x-z 
plane, every structure with q^ 1/2 is double-degenerate. 

We start with classification of the lattices exactly com- 
mensurate with the layered structure giving the set of 
commensurate fields (see also Ref. 0). The analysis of 
commensurability conditions can be done most conve- 
niently in scaled coordinates (y, z). In these coordinates 
the ground-state configuration corresponds to regular tri- 
angular lattice with period ao = •y/2$o7\/^7S^S^- It is 
convenient to consider orientation of the layered struc- 
ture with respect to this lattice rather than the other 
way round. The layered structure fits this lattice only if 
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it runs along one of the crystallographic directions, see 
Fig. njj. The direction (n, m) is defined by the lattice 
vector, e„^m, which can be expanded over the two basic 
lattice vectors, e(-„ — nei + me2. For nonequivalent 
directions n and m must be relatively prime numbers. 
In particular, two aligned configurations correspond to 
{n,m) = (1,0) and (1,1). Any such direction corre- 
sponds to set of matching fields, B(^n,m){N), which can 
be found by direct geometrical calculation^ 
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This result essentially rely on the London approx- 
imation, which implies a very strong inequality 
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Tin? 3> 1. Number of vortex- free layers 



per unit cell is given by — 1. The case N ~ 1 rep- 
resents a special situation when all the layers are filled 
with vortices and equivalent. It is interesting to note 
that even for dilute lattice one can have the Josephson 
vortices in every layer [N = 1) in the case of high-order 
commensurability (n, rn 3> 1). In ideal situation the lat- 
tice switches between different commensurate configura- 
tion via series of first-order phase transitions. Number 
of competing states rapidly increases with field decrease. 
In addition to giving general ground state, these lattices 
describe multiple metastable states with unique hierar- 
chical properties studied in Refs. 013 

Full analysis of the structure evolution requires energy 
consideration. The London model does not not describe 
layered superconductor at high fields. To obtain lattice 
structures in this region one has to consider more general 
Lawrence-Doniach model. The transition between the 
aligned lattices have been studied within this model by 
Ichioka^. However, our analysis shows that at many fields 
the true ground state is not given by aligned lattice. If we 
limit ourself only to aligned configurations, we reproduce 
results of Ref. at high anisotropies. 

Within the Lawrence-Doniach model, at fields ^ 
$o/(47rAAc) the lattice energy can be represented as 
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-u{N,q,h) 
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where A = Aa& and Ac are the components of the Lon- 
don penetration depth, h = 2ttjs^Bx/^o is the reduced 
magnetic field, and the reduced energy u{N, q, h) is given 
by 



(3) 



-1-1 - cos {(j)n+l - (f)n - hy)) 



To match the London limili^, we write u{N, q, h) in the 
form 

u{N,q,h)^^ln^ + lA32 + G{N,q,h) (4) 




FIG. 2: The field dependence of the energy function 
G{N, q, h) for the ground state and competing states. Each 
branch corresponds to the minimum of G{N, q, h) with re- 
spect to q at fixed h and A^. The curves are marked by value 
of A*'. Lattice configurations in scaled coordinates are shown 
at six marked fields. The left inset shows the field dependence 
of the shift parameter q near the structural transition of the 
dense lattice. The right inset illustrates competition between 
different configuration at smaller fields. The low-field local 
minima are marked by the indices corresponding to commen- 
surate configurations in the format {n,m),N. 



where the function G{N, q, h) defined by this equation 
approaches the London limi^, GL{r = iV^/i/(27r), g), for 
/i — > with 



GL{r,q) 



irr 1 exp (— 27rr/) -|- cos {2TTql) 1 

T"^^ 1 cosh (27rrZ) - cos {2nql) ~ 2 



This function depends only on lattice shape. Its absolute 
minimum corresponding to the triangular lattice is given 
by Gl(V3/2, 1/2) = -0.4022. 

We explore the evolution of the ground-state config- 
uration by direct numerical minimization of the energy 
© with respect to lattice parameter N and q defined 
in Fig.^ i.e., we computed reduced ground-state energy 
defined as G{h) = minN^q[G{N,q,h)]. The field depen- 
dence of the energy function G{N, q, h) is shown in Fig. 
|21for the ground state and competing states. Each curve 
corresponds to the minimum of G{N, q, h) with respect 
to q at fixed h and N and it is marked by the value of 
N. Every branch crossing corresponds to a first order 
phase transition between different commensurate states. 
Below, we show the first six lattice configurations which 
are realized with field decrease. At low fields we specify 
for each of these ground-state configurations correspond- 
ing indices (n, m) and period N in the format (n, m), A'^. 
At small fields the local minima of branches occur at 
commensurate fields corresponding to 

Let us review the lattice evolution with decreasing 
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magnetic field. At high fields an aligned dense Joseph- 
son vortex lattice is realized (structure a). This lattice 
becomes unstable at ft, w 1.33. Below this value, shear 
deformation develops corresponding to decrease of q be- 
low 1/2 (see structure b). The field dependence of q near 
the transition point is shown in the left inset of Fig. 
13 At /i « 0.99 this sheared dense lattice is replaced 
by the aligned lattice with the period N — 2 (structure 
c). At /i « 0.8 this structure transforms back into the 
misaligned lattice with period N = 1 (structure d). At 
smaller fields, many lattice configurations compete for 
the ground state, as one can see more clear in the right 
inset. As a consequence transitions become more and 
more frequent at smaller field. At several fields (e.g., at 
h «, 0.19, 0.137, 0.105 . . .) one or more lattice configura- 
tions have energies very close to the ground-state energy. 
We also note that there are several extended field ranges 
where in the ground state all layers homogeneously filled 
with vortices {N — 1) even in the region of dilute vor- 
tex lattice, e.g., for 0.115 < h < 0.17, 0.21 < h < 0.38. 
The layered structure favors such states. We also found 



that the layered structure does not favor aligned struc- 
tures with indices (n, m) — (1, 0). Such structures do not 
realize in ground state for 2 < N < 7. 

In conclusions, we explored ground states of the 
Josephson vortex lattice in layered superconductors. 
With decreasing field a dense lattice transforms into a 
dilute lattice via an intermediate sheared dense lattice 
state. After that the lattice goes through sequence of 
states with different orientations and c-axis periods sepa- 
rated by first-order phase transitions. At low fields many 
lattice configurations compete for ground state and phase 
transitions become very frequent. As the energy dif- 
ferences between different states induced by the layered 
structure become tiny at small fields, external factors, 
such as interactions with boundaries or with correlated 
disorder, may play role in selection of ground states in 
real samples. Frequently such external interactions favor 
configurations aligned with the layered structure. 
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